Magnetic Collapse of a Neutron Gas: No Magnetar Formation 
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A degenerate neutron gas in equilibrium with a background of electrons and protons in a magnetic 
field exerts its pressure anisotropically, having a smaller value perpendicular than along the magnetic 
field. For critical fields the magnetic pressure may produce the vanishing of the equatorial pressure 
of the neutron gas, and the outcome could be a transverse collapse of the star. This fixes a limit to 
the fields to be observable in stable pulsars as a function of their density. The final structure left 



over after the implosion might be a mixed phase of nucleons and meson (tt ± '°, 
strange star also likely) or a black string, but no magnetar at all. 



K ' ) condensate (a 



Gravitational collapse occurs in a body of mass M and 
radius R when its rest energy is of the same order of its 
gravitational energy, Mc 2 ~ GM 2 /R. We would like to 
argue that for a macroscopic magnetized body, e.g., com- 
posed by neutrons in an external field B, new physics 
appears and a sort of collapse occurs, when its internal 
energy density U is of the same order than its magnetic 
energy density M.B, where M. is the magnetization. This 
problem is interesting in the context of astrophysical and 
cosmological objects, as neutron stars. Extremely mag- 
netized neutron stars, or magnetars have been at- 
tracting the attention of reseachers in the last years . 
Their magnetic fields are estimated to be of order of 10 15 
G ||. For fields of this order of magnitude, there are 
values of the density for which the magnetic energy of 
the system becomes of the same order than its total en- 
ergy. For these physical conditions, superdense matter 
composed of neutral particles having a magnetic moment 
may undergo a transverse collapse since its pressure per- 
pendicular to B vanishes. This implosion is driven by the 
same mechanism described in Q for charged particles. 
As discussed below, the resulting object may be a hybrid 
or a strange star (depending on the macroscopic proper- 
ties of the imploding neutron star) but no any magne- 
tar. In the present paper (a preliminar report of which 
was presented in j^]) we will approach this problem by 
starting from a model of a relativistic degenerate neu- 
tron gas, pervaded by an extremely strong magnetic field, 
in equilibrium with a background of electrons and pro- 
tons. The last ones needed to prevent neutron beta decay 
n — > p + e + v through Pauli's principle. This configu- 
ration is maintained in approximate equilibrium through 
the equation ji n = fi p + fj, e + fip among their chemical po- 
tentials. However, as there is a flux of neutrinos escaping 
from the star its chemical potential would be taken as 
zero. 

The Dirac equation for a charged particle having an 
anomalous magnetic moment and placed in a strong mag- 



netic field was first derived by Pauli |7j . The free particle 
spectrum was obtained by Ternov et. al || . For neutral 
particles, it was presented first by Lostukov and Leven- 
tuiev, and reported by Klimenko We get for neutrons 
the eigenvalues, 
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where p%, p± are respectively the momentum compo- 
nents along and perpendicular to the external field B, 
q = 1.91Af„, where M n is the nuclear magneton, r\ — 
1,-1 are the 03 eigenvalues corresponding to the two 
orientations of the magnetic moment (parallel and an- 
tiparallel) with regard to the field B. The expression 
(0) shows manifestly the change of spherical to axial 
symmetry with regard to momentum components. In 
the non-relativistic limit and for strong fields, this is re- 
flected as small corrections to the spectrum as E(p, B) = 
p 2 /2m n + riqB{\ + p 2 ± /2m 2 l ) + q 2 B 2 /2m n , valid for any 
neutral Fermi gas with zero electric dipole moment and 
nonzero magnetic moment. 

The partition function Z which is obtained from the 
density matrix describing the model, leads to a thermo- 
dynamic potential Q = —T\nZ involving the contribu- 
tions from the species involved, which are considered to 
be in chemical equilibrium among themselves. 

Having an equation relating the chemical potentials, 
and demanding conservation of both baryonic number 
N n + N p = Nb and electric charge, N p + N e = 0, one 
may think to solve exactly the problem in terms of the 
external field as a parameter. In the present paper, how- 
ever, we focus our discussion on the properties of the 
equation of state.0 As it is known, the denser the Fermi 
gas, the better the ideal gas approximation jlCj ]. This 
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property is still valid in presence of external fields. We 
thus proceed to calculate the ideal gas thermodynamical 
quantities, and begin with the neutron thermodynamic 
potential. 

For the one- loop thermodynamical potential, f2„, we 
can write the expression 



f2„ — — f2o J]] [ 
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field at which the coupling energy of the anomalous mag- 
netic moment equals the rest energy is B cn = 1.57 x 10 20 . 
For protons B cp — 2.29 x 10 20 G, while for electrons it is of 
order B ~ e/rg ~ 10 _16 G, being ro the classical electron 
radius. However, these estimates must change if correc- 
tions depending on the field are made to the anomalous 
magnetic moment. At these critical fields vacuum be- 
comes unstable under particle-antiparticle pair produc- 
tion. Such decays are, however, suppressed in the dense 
gases by Pauli's principle. 

By defining x p = [i p /m p , y p = q p /m p , b = 2e/m 2 



s(x,r)y)], (2) q p = 2.79M„, then y p = 2.79e/2m 2 . We name also 



Above we have defined the variables x — i-i n /m n and 
y = qB/m n , and the functions 



f(x,w) = \Jx 2 - {l + i]y) 2 , 
s{x,rjy) = (vr/2 - sm _1 (l + T)y)/x, 
L(x,T)y) = ln(x + f(x,r)y))/(l + rjy). 



(3) 



For the density, defined as N n — —dCl n /dp, n , we can 
write 



N n = N [ 



f(x,rjy) 3 r]y(l + r)y)f(x,w) 
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while for the magnetization, given by M, 
-dQ n /dB, we get 



M n = -Mq 

•7=1,-1 



(1 - 2rjy)xf(x,Tiy) 



g(x p , B,n) = J x p — h(B,n) 2 , and 



h(B, n) = {y/bBn + 1 + rjy p B), (6) 
to have for the star's proton thermodynamical potential 
eBm 2 



n ±r/ 



x In 



x p +g(x p ,B,n) 



h{B,n) 
and for its density 

n ±77 

The magnetization is given by 

2 

M p = ^2^2[x p g{x p , B, n) - [h{B, nf + (r)y p 
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L(x,r)y) + —s(x,rjy)], (5) 



where N = M 3 /4tt 2 ~ 2.04 x 10 39 , Q = N M n ~ 
3.0 x 10 36 , and M = M 2 e^/4w 2 ~ 2.92 x 10 16 and one 
can write M n = M + {-q = -1) - M+(r) = +1), and 
obviously, M n > 0. 

If we include the anomalous magnetic moment for elec- 
trons, one can give a common formula for the spec- 
trum of electrons and protons in the external field B 



as |§: E = ^Jpl + (^j2eBn + ml p + r)q e , p B) 2 , where 
q e = ae/8nm e , q p ~ 2.79M„. For neutrons, the critical 



cent paper [[L5| on neutron gas in a magnetic field, containing 
expressions similar to ours for the spectra and densities of 
neutrons and protons (already reported in M), but different 
equations of state. The behavior of proton-neutron fraction 
and pion condensation are studied in detail in that paper. 



(bn/2VbBn+ l))]ln- 



+ \J x l~ g(x P ,B,n) 
h(B,n) ' 



(9) 



where the coefficents of these formulae are No = 
em p B/2ir 2 ~ 4.06 x 10 19 B, fl Q = N m p B ~ 6.1 x 10 16 B, 
and Mo = Nom p = Qo/B. 

For low magnetic fields n max is very large, and one can 
approximate the sum over Landau quantum numbers by 
an integral. This would lead to a neutron to proton ratio 
N n /N p > 8, similarly as the zero magnetic field case ||. 
The maximum occupied Landau quantum number n may 



be given as n Tl 



mjpBf - i/bB. 



For B <C B cp , so that y p B <C 1, and x p > 1, one can 
take approximately n max ~ (x 2 p — l)/bB, and for fields 
large enough n max = 0. As i p ~ i„, the proton density 
decreases with increasing B, favoring the inverse beta 
decay. In the extreme case of confinement to the Landau 
ground state we can estimate the bound N p = N e < 
2.5 x 10 9 BA^ /3 . For this limit it is roughly n n /n p > 
4 x 10- 10 7V 2 / 3 /5. For fields B ~ m p /q p and x p > 1, 
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n-max > 1, and thus large Landau numbers are again 
occupied. However, the dominant longitudinal pressure 
and magnetization comes from the neutron gas. 

The lost of rotational symmetry of the particle spec- 
trum determines an anisotropy in the thermodynamic 
properties manifesting in different equations of state for 
directions parallel and perpendicular to the external field, 
as is seen from the energy-momentum tensor in the con- 
stant magnetic field T^ v (ll]] , pj . Its spatial components 
Tn = T-22 = P±_ and T33 = P3, which contain the sum of 
the partial pressures of the several species involved, are 



Pl = -SI - BM 



Pz = -n, 



(10) 



where SI = ^T Sli is the total thermodynamical potential, 
i = n,p, e, and M = JT -Mi. By calling S = J2i &i an d 
N = Ni, where both the partial entropy and density 
are Si = dSli/dT, Ni = dSli/dfii, one can write the inter- 
nal energy density as U — — T44 = pN + TS + SI (we take 
T = in the present degenerate case), where the total 
quantities U, fxN, SI are of similar order. For positive 
magnetization, which is our case, the transverse pres- 
sure exerted by the magnetized particles is smaller than 
the longitudinal one in the amount BM. If we assume 
that the body is in equilibrium under the balance of neu- 
tron and gravitational pressures, being the latter of order 
Pgrav ~ GM 2 /R 4 , where R is the geometric average ra- 
dius of the star, the body stretch along the direction of 
the magnetic field. (This effect can be figured out from 
looking at the spectrum described by Eq.(Q), since the 
contribution from 77 = —1 terms are dominant. Then, 
if one approximates the second term inside the square 
root as [m n + (p 2 L /2m n — qB)] 2 , the term in parenthe- 
sis accounts for the transverse kinetic energy. This term 
decreases as the magnetic field increases. Equating it to 
zero, by taking p± ~ pp, where pp is the Fermi momen- 
tum, we obtain a functional relation between p n and B 
leading to the vanishing of the transverse kinetic energy, 
and in consequence of the transverse pressure. Observe 
that for x = 1.001, y~10 -6 . A more accurate result is 
obtained, however, from Tj_ = 0). This extreme case of 
vanishing of P± = —CI — BM. means that the transverse 
gravitational and Fermi gas pressures cannot compen- 
sate each other and an instability appears leading to a 
transverse collapse. We do not enter in the quantitative 
study of this problem, which would lead to a special sort 
of hybrid stars @H) or black strings (130]. In Fig- 
ure 1 we have drawn the equation P± n = in terms of 
the variables x,y. We observe that there is a contin- 
uous range of values of the chemical potential, starting 
from values x — 1.001 and magnetic field intensities, from 
y = 0.000001, for which the collapse takes place. The lat- 
ter value of y means fields on order of B ~ 10 14 G. To 
these ranges of x, y corresponds a continuous range of 
densities, from 10 _4 iVo (10 11 g/cm~ 3 ) onwards. Thus, 
although stable neutron stars start to occur for densi- 
ties N ~ 10 _2 iVo E^], the mechanism of decreasing the 
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FIG. 1. The instability condition: P±_ = 0. A neutron star 
having a configuration such that its dynamical stage would 
be represented by a point above the central curve in this plot 
would be unstable to transversal collapse, since P± < there. 



transverse pressure may act in the external regions of 
them, where densities of order 10 _4 iVo may occur. This 
would lead to transverse compression of the whole mass 
of the star, concentrating the magnetic lines of force and 



increasing B. For B ~ B c 



the magnetic coupling 



of quarks with B becomes of the order of their binding 
energy through the colour field. This might lead to a 
deconfinement phase transition leading to a quark (q)- 
star, a pressure-induced transition to uds-quark matter 
via ud-quark condensates, as discussed in Refs.( |l8|jL9| ] ) . 

Observations of GRB980827 from SGR 1900+14 ev- 
idenced the existence of a stable pulsation with period 
5.16 s [|. In analogy to the case for SGR 1806-20 
the observed spin-down rate of the pulse period, 
P = 1.1 x lO" 10 s s" 1 , led Kouveliotou et al. [g] to an- 
nouncing the discovery of a magnetar in the source SGR 
1900+14. These findings apparently lead to verify the 
Duncan and Thompson [jj magnetar model for SGRs. 
Kouveliotou et al. || claimed that spin down may be ex- 
plained by the emission of dipolar radiation from a NS 
endowed with a very strong magnetic field B ~ [2— 8]10 14 
G. Next we briefly review the standard theory of magne- 
tars. We shall show why they cannot survive after reach- 
ing the claimed superstrong magnetic fields, and then we 
present prospectives for a hybrid or strange star to ap- 
pear as a remnant of the quantum collapse. According to 
Duncan and Thompson [0 , neutron stars (NSs) with very 
high dipole magnetic field strength, Bp, ~ [10 14 — 10 15 ] 
G, may form when (classical) conditions for a helical dy- 
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namo action are efficiently met during the seconds fol- 
lowing the core-collapse in a supernova (SN) explosion 
E|. A newly-born NS may undergo vigorous convection 
during the first 30 s following its formation. If the NS 
spins (differentially) sufficiently fast (P ~ 1ms) the con- 
ditions are created for the a — ft dynamo action to be 
built. Collapse theory shows that some presupernova 
stellar cores could endow enough spin so as to rotate near 
their Keplerian equatorial velocity (the break-up spin) 
> ([§] 3 GAf/i? 3 ) 1/2 after core bounce. Under these 
conditions, fields as large as B ~ 10 17 (i£i) G may be 
generated as long as the differential rotation is dragged 
out by the growing magnetic stresses JlJ . For this process 
to efficiently operate the ratio between the spin rate (P) 
and the convection overturn time scale (r co „), the Rossby 
number R , should be < 1 (i?o ^> 1 should induce less 
effective mean-dynamos |ij). In this case, amplification 
of the magnetic field strength by these helical motions is 
not precluded, since the a 2 or a — Q dynamos may sur- 
vive depletion due to turbulent diffusion. An ordinary 
dipole B sat ~ [10 12 — 10 13 ] G may be built by incoher- 
ent superposition of small dipoles of characteristic size 
A ~ [| - l]km and B sat = (Airp) 1 / 2 \/r con ~ 10 16 G. At 
such fields, the huge rotational energy of a / > 1kHz NS 
is leaked out via magnetic braking, and an enormous en- 
ergy is injected into the SN remnant. This energy may 
power a plerion in the SN remnant. 

As shown above, at the end we are left with a NS 
with an extremely high field strength and a huge matter 
density p ~ [10 14 — 10 12 ] gem -3 . As illustrated in Figure 
1, those are the conditions for the quantum instability to 
start to dominate the dynamics of the young pulsar. At 
this stage, the magnetic pressure inwards may overpass 
the star's energy density at its equator, as defined in 
Eq. (|l0|) , and the collapse becomes unavoidable. As the 
collapse proceeds, higher and higher densities are reached 
till the point the supranuclear density may reverse the 
direction of implosion (hybrid star). From that moment, 
the sound wave generated at the core bounce builds itself 
into a shock wave traveling through the star at Vsw = 
c/ y/3 kms -1 . Although the magnetic field strength could 
be quite large as the collapse advances, the kinetic energy 
(E ~ 10 51 erg, the mean energy obtained in simulations 
of SN driven by the prompt shock fllfj) carried away by 
the shock wave may counterbalance it, and even surpase 
it, i.e., its ram pressure will equal the magnetic pressure 
at ta, the Alfven radius given by 



B 2 



8-Kpo 



2vr 2 



Gp. 



1/7 



Mm 2 



1/7 



(11) 



Then the magnetic field lines are pushed out, and fi- 
nally broken, from ta onwards, into the SN remnant sur- 
roundings as a violent explosion that dissipates a large 
part of the magnetic flux ~ B 2 r\) trapped inside 



the magnetar magnetosphere. This is analogous to the 
mechanism operating during a solar flare or a coronal 
mass-ejection, where the very high B in the "sun-spot" 
is drastically diminished. Although the process is quite 
fast, the large amount of matter ejected from the star 
at such large velocities drains out the dipole field of 
the remnant below the quantum electrodynamic limit of 
Bqed ~ 10 13 G. Further, since all the differential rota- 
tion has been dragged up to build up the superstrong 
magnetic field, then nothing else remains to make it to 
grow to its precollapse value. Thus no such ultra high B 
should reappear. We may be left with a submillisecond 
strange star JlJ] or a hybrid star (ll| with "canonical" 
field strength. 

We conclude that if a degenerate neutron gas is under 
the action of a strong magnetic field B < B c , for adequate 
values of its density its transverse pressure vanishes, the 
outcome being a transverse collapse. This phenomenon 
establishes a limit to the magnetic field expected to be 
observable in a neutron star, as a function of its den- 
sity, and suggests a possible end in the evolution of the 
highly magnetized neutron star as a mixed phase of nu- 
cleons and tt^ 1 ' , k^' meson condensate or a black string. 
Thus this result implies that most of the observed pul- 
sars should have some either strange matter or meson 
condensate in their interiors. 
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